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ABSTRACT 

We study the effects of a dark energy component with equation of state p — wp with constant w > — 1 on the 
formation of Cold Dark Matter (CDM) haloes. We find two main effects: first, haloes form earlier as w increases, 
and second, the amplitude of the dark-matter power spectrum gets reduced in order to remain compatible with the 
large scale Cosmic Microwave Background (CMB) anisotropies. These effects counteract. Using recipes derived 
from numerical simulations, we show that haloes are expected to be up to ~ 50% more concentrated in CDM 
models with quintessence compared to ACDM models, the maximum increase being reached for w ~ —0.6. For 
larger w, the amplitude of the power spectrum decreases rapidly and makes expected halo concentrations drop. 
Halo detections through weak gravitational lensing are highly sensitive to halo concentrations. We show that weak- 
lensing halo counts with the aperture-mass technique can increase by a factor of ~ 2 as w is increased from —1 to 
—0.6, offering a new method for constraining the nature of dark energy. 



1. Introduction 
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5 In recent years, cosmology has seen increasing observational ev- 
£^ idence for an accelerating phase of the cosmic expansion, most 
notably through the observations of distant type la supernovae 
C"| (Perlmutter et al. 1999, Riess et al. 1998). This astonishing evi- 
O dence motivated renewed interest in the properties of the energy 
density ascribed to the "vacuum". A vacuum energy compo- 
nent should account for both the accelerating expansion and for 
the residual ~ 70% of the energy density required for recon- 
^ ciling the geometrical flatness required by Cosmic Microwave 
Background (CMB) observations (De Bernardis et al. 2002, Lee 
et al. 2001, Halverson et al. 2002) with the evidence of a low- 
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density universe with Ho ^ 0.3 (Percival et al. 2001). 
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While one of the historical candidates for such an energy den- 
sity is the cosmological constant, introduced as a simple geomet- 
rical term in Einstein's equations, it is well known that it leads 
to serious and unsolved theoretical problems. The exceedingly 
low value of the vacuum energy density today, compared to that 
allowed by the most plausible theories of the early stage of the 
Universe, motivated the introduction of a more general concept 
now widely known as "dark energy". 

Preceding the evidence for cosmic acceleration, a generali- 
sation of the cosmological constant by means of a scalar field, 
now known as "quintessence", was proposed (Wetterich 1988, 
Ratra & Peebles 1988). In this class of models, the dark energy 
is supposed to reside mostly in the potential energy of the field, 
which interacts only gravitationally with ordinary matter. The 
evolution is described by the ordinary Klein-Gordon equation. 
If the potential is flat enough, or if the motion of the field along 
its trajectory is sufficiently slow, a cosmological constant-like 
behaviour can be mimicked by the scaling of the energy density 
of the quintessence field. 

For general forms of the scalar field potential, there exist at- 
tractor trajectories for the evolution of the background expecta- 
tion value of the field. These trajectories are known as "track- 
ing" (Steinhardt, Wang & Zlatev 1999) and "scaling" (Liddle & 
Scherrer 1999) solutions. They have been shown to alleviate, 
at least at a classical level, the fine-tuning required in the early 
Universe, when the typical energy scales were presumably com- 



parable to the Planck scale, to generate a vanishing relic vacuum 
energy as it is observed today, 120 orders of magnitude smaller. 
However, these scenarios are not able to solve the coincidence 
problem, i.e. why we are living in the epoch in which dark 
energy and matter have roughly the same energy density. De- 
spite some attempts at addressing this issue (Tocchini-Valentini 
& Amendola 2002, Chiba 2001, Armendariz-Picon, Mukhanov 
& Steinhardt 2001, Dodelson, Kaplinghat & Stewart 2000), it 
remains one of the greatest puzzles of modern cosmology. 

For constraining the nature of the dark energy, an important 
step would be accomplished if parameters could be constrained 
which capture its most essential features. In particular, if the 
dark energy is modelled as a quintessence field in the tracking 
regime, the simplest description of its properties will require the 
use of only two parameters, i.e. its present energy density and 
the ratio between pressure and energy density in its equation of 
state. Generally, this ratio depends on time, as implied by the 
evolution according to the equation of motion. However, it can 
be shown that, in most simple models of quintessence involv- 
ing an inverse power-law potential, the effect of a time variation 
of the equation of state can be neglected at low redshifts, when 
the field has settled on its tracking trajectory. In this case, the 
equation of state is simply related to the power with which the 
potential depends on the field itself. This simplification allows 
constructing a scheme for describing the dark energy behaviour 
at redshifts where interesting cosmological effects arise, such 
as the effect on the magnitude-redshift relation of type la su- 
pernovae (Perlmutter 1999, Riess 1998) and the effect on grav- 
itational lensing of distant galaxies and quasars (Futamase & 
Yoshida2001). 

Even though the dark energy dynamics has a geometrical ef- 
fect on acoustic features of the CMB anisotropy (Baccigalupi 
et al. 2002; Doran, Lilley & Wetterich 2002; Corasaniti & 
Copeland 2002), it is now commonly accepted that the most 
interesting properties of a dark energy component have to be 
probed by looking at processes occurring at low redshifts when 
it starts dominating the cosmic expansion. 

Thus, one of the most powerful probes of the quintessence 
field results from its effects on cosmic structure formation, most 
notably on the background cosmology and the evolution of indi- 



vidual collapsing overdensities. First, a dark energy component 
affects the matter density of the background in which haloes 
form. Second, the amplitude of matter perturbations is sensi- 
tive to the presence of a dynamical vacuum energy, through the 
normalisation of the matter power spectrum to the large, un- 
processed, scales probed by the large-scale CMB anisotropies. 
Third, changes in the background matter density induced by a 
dark energy component can seriously affect characteristic prop- 
erties of collapsing structures. As shown by Lokas & Hoffmann 
(2002), a substantial quintessence component changes the char- 
acteristic density of a forming dark matter halo. 

In this paper, we study how quintessence affects the con- 
centration of dark-matter haloes, and resulting changes in their 
weak-lensing efficiency. Weak gravitational lensing provides a 
powerful tool for mapping the large-scale mass distribution (see 
Mellier 1999a and Bartelmann & Schneider 2001 for reviews), 
and the potential impact of dark energy on the weak lensing con- 
vergence power spectrum has already been recognised (see Mel- 
lier 1999b; Huterer 2002). We show in this paper that an energy 
density component with negative pressure affects weak lensing 
by dark-matter haloes not only through changes in the global 
properties of the Universe, but also by modifying their internal 
density concentration. The main idea is that dark energy affects 
structure growth and thus the time of halo formation. Since halo 
concentrations reflect the density of the Universe at their forma- 
tion epoch, this affects halo mass distributions, and weak lensing 
provides methods for quantifying respective changes. 

The paper is organised as follows. In Sect. 2, we describe the 
main effects of the dark energy on the cosmological growth fac- 
tor, volume elements, and the normalisation of the dark-matter 
power spectrum. In Sect. 3 we compute the impact on halo con- 
centration. In Sect. 4 we predict resulting effects on the weak- 
lensing aperture mass statistics, and Sect. 5 contains our conclu- 
sions. 



2. Cosmological implications of quintessence 

MODELS 

We model the dark energy as a spatially homogeneous compo- 
nent, labelled Q, with constant equation of state parameterised 
by the ratio between the pressure pq and the energy density Pq, 
w — pq/pq. We neglect a possible time variation of w , as well as 
any effects possibly due to spatial inhomogeneities of the dark 
energy. Indeed, at least in most models proposed so far, the 
relevant cosmological effects of the dark energy compared to 
a cosmological constant are mainly related to its effective equa- 
tion of state at redshifts when it is relevant for cosmic expansion, 
say z < 5, as we already noted in the introduction. Neglecting 
inhomogeneities of the dark energy is justified in the present 
context since they are likely to show relativistic behaviour on 
sub-horizon cosmological scales; indeed, the effective mass of 
the vacuum component, of the order of the critical density today, 
is extremely light compared to any other known massive parti- 
cle, so that quintessence clustering occurs only on scales larger 
than or equal to the horizon size (Ma et al. 1999). It can also be 
shown formally that a minimally coupled quintessence field has 
a relativistic effective sound speed (Hu 1998), so that its fluctu- 
ations are damped out on the scales in which we are interested 
here. 

Assuming the equation of state pq = wpQ, the adiabatic equa- 
tion implies 
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Friedmann's equation can then be written, 



Pq = Pq.oo 
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where Q.q is the quintessence density parameter at the present 
epoch, i.e. Pq.o divided by the critical energy density today. The 
density parameter for non-relativistic matter is i2o today. Ob- 
viously, for w = — 1, Friedmann's equation for a cosmological 
constant is retained. Unless stated otherwise, we will assume 
in this paper that the curvature of spatial hypersurfaces is zero, 
thus £2q + £2o = 1, and the curvature term in (Q) vanishes. 

The quintessence term in (^) has two immediate conse- 
quences relevant for our purposes. First, the way how density 
inhomogeneities grow is modified, and second, the cosmic vol- 
ume per unit redshift changes. 

In linear theory, the density contrast 8 of matter perturbations 
grows according to 
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(Heath 1977; see also chapter 15 of Peacock 1999). The right- 
hand side of (^|) is proportional to the so-called growth factor 
D + (a), which is commonly normalised either to unity at a = 1, 
or such that it rises proportional to a for a <C 1. We plot in 
Fig. [j] the growth factor as a function of redshift z, normalised 
to unity today and divided by a. The normalisation ensures that 
D + /a goes to unity for z — > 0. In an Einstein-de Sitter model, 
D + = a, thus the curves show how much faster structures grow 
in the assumed model universes compared to an Einstein-de Sit- 
ter model. Curves are shown for a variety of cosmological mod- 
els. All of them have Q.q —0.3 and either £2q = 1 — Q.q or, for 
comparison, Q.q = 0. 




FIG. 1 . — Growth factor D+ (z) as a function of redshift for five different 
cosmological models as indicated. The growth factor is normalised to 
unity at the present epoch, and divided by the scale factor to emphasise 
the differences between the models. With increasing w, the growth 
factor increases towards its value for the open model with Q.q = 0. 

The solid curve in Fig. [l] shows the growth factor for a model 
universe with a cosmological constant (w = — 1). Going back 
in time, D + /a rises from unity to ~ 1.3 and turns flat near 



3. This means that structures start forming more 



where Pq.o is the quintessence energy density today, and a is 
the cosmic scale factor normalised to unity at the present epoch. 



quickly in this model compared to an Einstein-de Sitter model, 
but the growth slows down considerably at redshifts smaller than 
z ~ 3. For the low-density open model without quintessence (la- 
belled £2q = 0), D + /a keeps rising as z increases. At redshift 
5, for instance, the amplitude of structures is 1.4 times higher 



than in the cosmological-constant model. Increasing w inter- 
polates between the open and the flat model with cosmological 
constant. Thus, keeping ilo and Q.q fixed, structures form ear- 
lier for larger values of w, approaching the growth behaviour for 
low-density open models without quintessence or cosmological 
constant. 

A similar interpolation is seen in the behaviour of the cosmic 
volume per unit redshift. Figure shows 



V(z) = / dz'4KD 2 (z') 
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where D(z) is the angular-diameter distance between redshifts 
and z, thus V(z) is the proper volume of a sphere of "radius" 
z around the observer. Similar to the behaviour of the growth 
factor, dark-energy models interpolate between the two limit- 
ing curves, where the model with cosmological constant has 
the largest and the model without cosmological constant has the 
smallest volume. 
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FIG. 2. — Cosmic volume in units of the Hubble volume between red- 
shift zero and z for six different cosmological models as indicated. With 
increasing w, the cosmic volume decreases towards its value for the 
open model with Q.q = 0. 

Here and below, we keep the cosmological parameters fixed 
at Q.Q = 0.3, Q.q = 0.7. The Hubble constant is Hq = 
lOO/ikms-'Mpc- 1 with h = 0.7. 

In order to describe the formation of dark-matter haloes, we 
need to specify the power spectrum of dark-matter fluctuations. 
We choose the cold dark matter (CDM) model, whose transfer 
function was given by Bardeen et al. (1986), and we set the index 
of the primordial power spectrum to the Harrison-Zel'dovich 
value of n = I. Then, the power spectrum has two free param- 
eters, the shape parameter F which locates its maximum, and 
the amplitude. For the shape parameter, we assume F = H^h 
as suggested by theory and in agreement with observations of 
the galaxy power spectrum. For our model for dark energy, the 
transfer function by Bardeen et al. is applicable because dark 
energy does not cluster on the relevant scales, and the shape pa- 
rameter r is set by the scale of matter-radiation equality, which 
is unaffected by w (cf. Wang & Steinhardt 1998). 

The amplitude of the power spectrum needs to be chosen 
such that certain observations can be reproduced. It is an im- 
portant constraint that the abundance of rich galaxy clusters in 
our cosmic neighbourhood be reproduced. Since the galaxy- 
cluster mass function falls very steeply at the high-mass end, 
small changes in the amplitude of the power spectrum lead to 



large changes in the cluster number density, thus the amplitude 
is in principle well constrained by the cluster abundance. How- 
ever, for that normalisation procedure, rich galaxy clusters are 
identified by their X-ray emission, thus the reliability of the nor- 
malisation depends on how well models can describe the X-ray 
properties of the intracluster gas. 

We saw before that structure starts forming earlier in 
quintessence models with w > — 1 compared to models with cos- 
mological constant. Galaxy clusters forming earlier are hotter 
because of the higher mean density of their surroundings. Re- 
producing the current number density of X-ray selected clusters 
thus requires a power-spectrum amplitude which decreases with 
increasing w (Wang & Steinhardt 1998). Figure| shows their re- 
sults, expressing the power-spectrum amplitude in terms of the 
rms fluctuation amplitude Og on a physical scale of 8/z _1 Mpc. 
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FIG. 3. — The normalisation of the power spectrum, expressed in terms 
of Og, is shown as a function of w for two different normalisation meth- 
ods. The solid line shows the cluster-abundance normalisation derived 
by Wang & Steinhardt, the dotted line shows as constrained by the 
COBE-DMR data. While rjg is almost constant in the latter case for 
w < —0.6, it drops rapidly for larger w because of the increasingly 
strong integrated Sachs-Wolfe effect. The uncertainty of cjs determined 
from COBE data is approximately 7% (Bunn & White 1997). 

The power spectrum of density perturbations can be nor- 
malised by the CMB anisotropies on large angular scales, cor- 
responding to physical scales which were larger than the hori- 
zon at decoupling. The only available data on those scales 
are those of the Differential Microwave Radiometer (DMR) 
experiment on board the COsmic Background Explorer satel- 
lite (COBE, Bennett et al. 1996), but much improved data 
are being taken by the M icrowave Anisotropy Probe (MAP, 
http: //map.gsfc.nasa.gov/), and will be taken by the Planck satel 



lite ( ^ittp://astro.estec.esa'nl/SA-generaFProjects/Planck| ). Sub- 
degree CMB anisotropies have been measured by several exper- 
iments (De Bernardis et al. 2001, Lee et al. 2001, Halverson et 
al. 2002). However, on these angular scales, the main effect of 
increasing w above — 1 is to rigidly shift acoustic peaks toward 
larger angular scales (see Baccigalupi et al. 2002 and references 
therein), having almost no effect on the overall normalisation of 
the spectrum. 

Normalising the power spectrum according to the COBE mea- 
surements fixes the amplitude of the power spectrum on its 
large-scale end. As w increases at fixed £2o and £2q, dark en- 
ergy dominates cosmic expansion earlier compared to a model 
with cosmological constant having the same energy density to- 
day. This has the effect of enhancing the dynamics of the grav- 
itational potential due to the change in the cosmic equation of 



state, thus increasing the Integrated Sachs-Wolfe (ISW) effect 
on COBE scales. While this trend is gentle for w < —0.6, it 
steepens for larger w, as shown by the dotted line in Fig. ||. 

This non-linear behaviour of ag reflects the sensitivity of the 
ISW effect (and thus the normalisation) to the redshift zq at 
which dark energy starts dominating over matter. It is easy 
to show that this is determined by the power law 1 + zq = 

(Q.q/Q.q)- 1 ^ w . We recall that, although we consider values of w 
up to —0.4 for the sake of generality, values larger than —0.6 are 
not interesting in the framework of dark-energy models. Indeed, 
w larger than —0.6 is insufficient to provide cosmic acceleration 
if Q.q = 0.7 as we assume in this work. In addition, the sharp de- 
crease in Fig. |]due to the ISW leads to a decrease of the acoustic 
peaks in the CMB power spectrum below the level observed by 
experiments operating on sub-degree angular scales. 

In order to avoid the uncertainties in Og due to the uncertain- 
ties in modelling the X-ray cluster population, we choose the 
COBE normalisation for our study. We adopt the normalisation 
method by Bunn & White (1997), which has a 7% accuracy. 
This procedure exploits a maximum likelihood approach for re- 
producing the measured CMB anisotropy power once the sky 
regions affected by the Galactic signal have been cut out. The 
CMB anisotropy can be expressed as a line-of-sight integral of 
the perturbation power spectrum weighted with suitable geomet- 
rical functions, implemented in the CMBFAST code (Seljak & 
Zaldarriaga 1996). For taking the dark energy component into 
account, we use a modified version of CMBFAST (see Bacci- 
galupi, Matarrese & Perrotta 2000). Since halo properties are 
determined by the small-scale end of power spectrum while the 
COBE normalisation fixes its large-scale end, quantitative re- 
sults will sensitively depend on the index n of the power spec- 
trum. 

An entirely different and perhaps more direct way of nor- 
malising the power spectrum has become feasible in the recent 
past. Large-scale density fluctuations differentially deflect light 
on its way from distant sources to us. The gravitational tidal 
field of the matter inhomogeneities coherently distorts the im- 
ages of faint background galaxies. Albeit weak, this cosmic 
shear effect has recently been measured successfully by several 
groups, whose results agree impressively although different tele- 
scopes, observational parameters and data-reduction techniques 
were used (Van Waerbeke et al. 2000; Bacon, Refregier & El- 
lis 2000; Kaiser, Wilson & Luppino 2000; Wittman et al. 2000; 
Maoli et al. 2001; Van Waerbeke et al. 2001). Since gravitational 
lensing depends only on the matter distribution and not on its 
composition or physical state, cosmic shear should provide one 
of the cleanest ways for constraining the power-spectrum ampli- 
tude. 

Cosmic-shear measurements have been shown to agree very 
well with the expectations in a ACDM universe. Fixing the 
shape parameter to Y = 0.2 1 , the measurements require Q.q < 0.4 
and Og > 0.7 at 95% confidence (Van Waerbeke et al. 2001). 
Note, however, that these constraints depend on the nonlinear 
evolution of the power spectrum which has some uncertainties. 

This result is almost completely insensitive to w in 
quintessence models. In Fig. f|, we show the rms cosmic shear 
in apertures of radius 9 as a function of 9 for four models 
differing by w, as indicated. The power spectrum was nor- 
malised to reproduce the COBE measurements, and its nonlinear 
evolution was approximated using the fitting formulae by Pea- 
cock & Dodds (1996). Following Van Waerbeke et al. (2001), 
we adopted the source redshift distribution shown in (||) below 
with parameters zq = 0.8 and p = 1.5. The curves are much 
closer than the typical uncertainty of cosmic-shear measure- 
ments, showing that the constraint on ag derived from cosmic 
shear is independent of w. The main reason is that cosmic shear 
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FIG. 4. — The squared cosmic shear { | Y 2 1 ) (©) ^ averaged in apertures of 
radius 0, is shown as a function of 9 for five different cosmological 
models as indicated. A source redshift distribution with mean redshift 
~ 0.9 was assumed for the plot. Changing w has very little effect on the 
curves, implying that the constraints on Og derived from cosmic-shear 
measurements are insensitive to w. 



is most sensitive to structures below redshift ~ 0.5 where the 
differences between the models are small. The figure shows that 
the COBE normalisation is not in conflict with the measured 
cosmic shear for all interesting values of w. 

3. Halo concentrations in dark energy models 

Numerical simulations of cosmic structure formation show con- 
sistently that the density profiles of dark-matter haloes can be 
described by a two-parameter family of models. Far outside a 
scale radius r s , the profiles fall off proportional to r~ 3 , while 
they are cuspy but considerably flatter well within r s . The exact 
inner profile slope is under debate. The second free parameter 
besides the scale radius is a characteristic density scale p s . We 
adopt the density profile suggested by Navarro et al. (1995), 



The two parameters r s and p s are not independent. Haloes are 
commonly parameterised by their virial mass M. For a given 
cosmology, this also defines their virial radius R. Numerical 
simulations show that the scale radius depends on the mass such 
that the halo concentration c =R/r s is a characteristic function 
of mass. Given M, R and c, the density scale p s is fixed. For 
definiteness, we parameterise halo masses consistently by M200, 
i.e. masses enclosed in spheres with radius R200 in which the 
average density is 200 times the critical density. 

Numerically simulated haloes turn out to be the more concen- 
trated the earlier they form. This is interpreted assuming that the 
central density of a halo reflects the mean cosmic density at the 
time when the halo formed. This implies that haloes forming 
earlier are expected to be more concentrated. 

Several algorithms based on this assumption have been sug- 
gested for describing the concentration of dark-matter haloes. 
Originally, Navarro et al. (1997) devised the following approach. 
A halo of mass M is first assigned a collapse redshift z co n de- 
fined as the redshift at which half of the final halo mass is con- 
tained in progenitors more massive than a fraction /nfw of the 
final mass. Then, the density scale of the halo is assumed to be 



some factor C times the mean cosmic density at the collapse red- 
shift. They recommended setting /nfw = 0.01 and C — 3 x 10 3 
because their numerically determined halo concentrations were 
well fit assuming these values. 

Bullock et al. (2001) noticed that halo concentrations change 
more rapidly with halo redshift than the approach by Navarro et 
al. (1997) predicts. They suggested a somewhat simpler algo- 
rithm. Haloes are assigned a collapse redshift defined such that 
the non-linear mass scale at that redshift is a fraction /b of the 
final halo mass. The halo concentration is then assumed to be a 
factor K times the ratio of the scale factors at the redshift when 
the halo is identified and at the collapse redshift. Comparing 
with numerical simulations, they found /b = 0.01 and K = 4. 

Yet another algorithm was suggested by Eke et al. (2001). 
They assigned the collapse redshift to a halo of mass M by re- 
quiring that the suitably defined amplitude of the linearly evolv- 
ing power spectrum at the mass scale M equals a constant C^ s . 
Numerical results are well represented setting Cens = 28. 
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FIG. 5. — Halo concentration parameters c as functions of halo mass 
M2oo- Results obtained from three different prescriptions for calcu- 
lating the concentration are shown; these are the prescriptions from 
Navarro, Frenk & White (NFW), Bullock et al. (B) and Eke, Navarro 
& Steinmetz (ENS). Curves are shown for the ACDM model (w = — 1) 
and for a quintessence model with w = — 1/2. Irrespective of the pre- 
scription, the concentration parameters are higher in the quintessence 
model than in the ACDM model for all halo masses. 

Since we consistently use /?200 and M200 for parameterising 
haloes, we need to convert masses and concentration parame- 
ters from the slightly different definitions introduced by Bullock 
et al. and Eke et al. In particular, this requires us to iteratively 
compute the concentration parameter according to our defini- 
tion. Moreover, the algorithms implicitly use the mean overden- 
sity of virialised haloes, A v , and the linear overdensity of col- 
lapsed haloes, 8 C . These parameters depend on cosmology and 
on the dark energy parameter w. We compute them using the 
formulae given in Lokas & Hoffmann (2002). 

Although halo concentrations produced by these different al- 
gorithms differ in detail, they have in common that haloes form- 
ing earlier are more concentrated. Taken together with the ear- 
lier result that haloes form earlier in dark energy models with 
w > — 1 than in ACDM models, this implies that haloes are ex- 
pected to be more concentrated in models with w > — 1. Fig- 
ure |] illustrates this. Halo concentrations computed with the 
three different algorithms are plotted as functions of halo mass 
for redshift zero. Since less massive haloes start forming earlier 
than more massive ones in hierarchical models like CDM, the 
concentration decreases with halo mass. Curves are shown for 



models with w = — 1 and w = — 1/2, keeping all other parame- 
ters fixed. AtM~5x 1O 13 /i _1 M , for instance, halo concentra- 
tions are approximately 50% higher in the dark energy compared 
to the ACDM model. 

Although concentrations of haloes identified at z > tend to 
be smaller than at redshift zero, the trend of increasing concen- 
trations with increasing w remains also at higher redshifts. 
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FIG. 6. — Halo concentrations obtained from the algorithm described 
by Bullock et al. as functions of halo mass M200 for five different cos- 
mological models as indicated. With increasing w, the concentration 
increases until w > —0.6, and drops rapidly as w increases further. 

Figure || shows halo concentrations as functions of mass de- 
rived from the algorithm suggested by Bullock et al. for five 
different choices of w, keeping all other parameters fixed. As 
w increases away from —1, concentrations first increase for all 
halo masses shown. A maximum is reached for w ~ —0.6. As 
w is raised further, concentrations decrease rapidly. This is an 
effect of the power-spectrum normalisation. As w increases, Os 
must decrease because otherwise the normalisation constraints 
would be violated; either, there would be too many hot galaxy 
clusters, or the secondary CMB anisotropies caused by the in- 
tegrated Sachs-Wolfe effect would exceed the COBE measure- 
ments. As Gg is lowered, haloes form later, thus counter-acting 
the earlier increase of the growth factor in quintessence mod- 
els. Adopting the COBE normalisation, the maximum effect is 
achieved just before the integrated Sachs-Wolfe effect sets in 
strongly at w > —0.6, cf. Fig. 

4. Weak lensing and halo counts 

The impact of dark energy on halo concentrations should cause 
an impact on many observable quantities. As an illustration, we 
will now describe how the number of haloes detectable through 
their weak gravitational lensing effect change with w. 

The gravitational tidal field of individual sufficiently massive 
haloes imprints coherent distortions on the images of faint back- 
ground galaxies in their neighbourhood. Haloes can thus be 
detected searching for their characteristic signature on the ap- 
pearance of the background galaxy population. A sensitive and 
convenient technique for halo detection, the aperture mass tech- 
nique, has been suggested by Schneider (1996). 

Consider a circular aperture of angular radius 9. The aperture 
mass is defined as a weighted integral of the lensing convergence 
K across the aperture, 

Map(6) = |d 2 #K(^M|£|). (6) 



The convergence is the surface mass density scaled by its critical 
value for strong lensing, 



K=— E ° s 
K E cr ' er 4kG D d D ds 



(7) 



where Aj.ds.s are the angular-diameter distance from the ob- 
server to the lens, the lens to the source, and the observer to 
the source, respectively. 

Through E C r, the aperture mass (Q) depends on the source red- 
shift. We compute mean aperture masses by averaging M ap over 
the normalised source-redshift distribution 



p(Zs) = 



P 



^r(3/p) 



exp 



(8) 



with zo = 1 and p = 1.5, implying a mean source redshift of 
- 1.5 (cf. Smail et al. 1995; Cohen et al. 2000). 
If the weight function t/(|r}|) is compensated, 



d 2 iW(|t}|) = 0. 



(9) 



the aperture mass M ap can be written as a (differently) weighted 
integral across the aperture of the tangential component of the 
shear with respect to the aperture centre. Thus, M ap is a directly 
observable quantity. 

Halo detection can then proceed as follows. An aperture of 
given radius is shifted across a wide and sufficiently deeply ob- 
served field. At all aperture positions, the aperture mass is deter- 
mined. Potential haloes are located where M ap exceeds a certain 
threshold. 

The signal-to-noise ratio of an aperture-mass measurement is 
given by S = M^o^ with the dispersion 



o M (6) =0.016 



30 arcmin 2 



-1/2 



\Q.2J \ 1' 



(10) 



where n g is the number density of faint background galaxies and 
G E is the variance of their intrinsic ellipticity distribution. We 
assume « g = 30 and take into account that only such sources 
from the distribution (^|) contribute to the signal whose redshift 
is larger than that of the lensing halo. 

Convolving the projected NFW density profile (cf. Bartel- 
mann 1996) with the weight function U the aperture mass 
of NFW haloes as a function of halo mass, M ap (M,0) is easily 
computed. Assuming further a halo mass function M, a io(yW,zX 
we can calculate the number of haloes per unit mass and redshift 
whose aperture mass is sufficiently high for the signal-to-noise 
ratio S to exceed a given threshold 5 m i n - We choose the mass 
function suggested by Sheth & Tormen (1999), which is a vari- 
ant of the Press-Schechter (1974) mass function which well re- 
produces the mass function found in numerical simulations. We 
take into account that our definition of mass differs slightly from 
Sheth & Tormen's in that we use the mass enclosed by a sphere 
in which the mean density is 200 times the critical rather than 
the mean density. 

In calculating the number density of haloes of mass M at 
redshift z which produce a significant weak-lensing signal, 
N] ens (M,z), we have to take into account that a signal-to-noise 
threshold 5 m j n for the weak-lensing signal does not correspond 
to an equally sharp threshold in halo mass because of the scat- 
ter in the aperture mass which is caused by the shot noise from 
the discrete background galaxy positions and their intrinsic el- 
lipticity distribution. A halo of mass M has a certain probability 



p(M ap |M) to produce an aperture mass M ap , which we model as 
Gaussian, 



/?(M ap |M)ocexp 



[M ap -M ap (M)] 
2a 2 , 



(11) 



where Om is given by ©and M ap (M) is the theoretical expec- 
tation for the aperture mass of a halo with mass M. The proba- 
bility for a halo of mass M to have an aperture mass above the 
signal-to-noise threshold is then 



P(S >5 min |M) = -erfc 



S{M) 



V2 



(12) 



where erfc(x) is the complementary error function and S (M) = 
M ap (M)oM . Thus, the number density of significantly lensing 



haloes is 



Aliens (M,z) =P(S > S min \M)N ha i (M,z) 



(13) 



We illustrate A^i ens (M, z) in Fig. ^| Here and below, we set 5 n 
5. 
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FIG. 7. — Contours in the mass-redshift plane showing the number den- 
sity of such haloes which are capable of producing a significant weak- 
lensing signal with the aperture-mass technique. The contour levels are 
drawn at a halo density of io~ 9 - 5 '~ 8 - 5 hM^ . The different line types 
indicate results for different cosmological models, as indicated. With 
increasing w < —0.6, contours widen, while they shrink as w increases 
further. 

The figure shows two sets of contours, inner and outer, whose 
levels are 10~ 8 5 and 10~ 9 5 hM- , respectively. Each set has 
four contours for different values of w, as indicated in the figure. 
All other cosmological parameters are kept fixed, i.e. Ho = 0.3 
and £2q = 0.7. The solid contour is for the ACDM model, for 
which w = — 1 . The dotted and short-dashed contours show that 
the region in the mass-redshift plane occupied by significantly 
lensing haloes widens as w increases from —1 to —0.6. If w is 
increased further, this region shrinks considerably, as the long- 
dashed contour shows. This illustrates the competition between 



the two effects outline above: haloes grow earlier and are thus 
more concentrated in dark-energy models with w > — 1, but the 
integrated Sachs-Wolfe effect reduces the power-spectrum nor- 
malisation required by the COBE-DMR data. The maximum 
extent of the contours in the mass-redshift plane is reached just 
before the strong decrease in Os near w ~ —0.6 illustrated in 
Fig- 1 
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FIG. 8. — The redshift distribution of weak-lensing haloes is shown for 
four different cosmological models, as indicated. As w changes, the 
peak remains near z ~ 0.4, but the amplitude increases until w ~ —0.6 
and drops for further increasing w. 

Figure || shows the redshift distribution of weak-lensing 
haloes for the same four dark energy models used for Fig. 0. The 
curves in Fig. |8] are thus integrals over mass of the distributions 
in Fig. 0, 

= /°°dMMens(M, Z ) . (14) 

az Jo 

As w increases from w = — 1 to w = —0.6, the amplitude of 
the redshift distribution rises by ~ 50%, and the redshift distri- 
bution extends towards higher redshift. At w = —0.4, the am- 
plitude drops to the level of the ACDM model (w — — 1), but 
the distribution is somewhat wider. This reflects the fact that 
the number of haloes is reduced compared to the models with 
somewhat lower w because the power-spectrum normalisation 
is lower, but the haloes are more concentrated compared to the 
ACDM model because of their earlier formation, allowing them 
to be significant lenses at higher redshifts. 

Finally, we show in Fig. ^ the total expected number of weak- 
lensing haloes per square degree as a function of w. The peak is 
reached with ~ 17 haloes per square degree at w — —0.6, which 
is almost a factor of two higher than for ACDM models. Note 
also that the increase is roughly linear with w up to the max- 
imum. For larger values of w, the halo number drops steeply. 
However, as already stressed, these cases are not interesting for 
cosmology since they do not produce cosmic acceleration as re- 
quired by observations. 

5. Summary and Conclusions 

We investigated the expected properties of dark-matter haloes in 
dark energy cosmologies. For our purposes, the essential fea- 
tures of such models are captured describing the dark energy as 
a density component with negative pressure pq — wPq, where 
w > — 1 is a constant. The dark energy density Pq is determined 
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FIG. 9. — Total number of significant weak-lensing haloes per square 
degree as a function of w. The curve reaches a peak near w ~ —0.6, 
where the halo number is approximately twice that for a ACDM model. 
For w >= —0.6, the halo number drops steeply. 



by its present value in units of the critical density, £Iq. In agree- 
ment with results from observations of the CMB, we focus on 
models which are spatially flat, Q.q + Q.q — 1 , have a matter den- 
sity parameter Q.q = 0.3, and a Hubble constant of h = 0.7. 

The modified background dynamics in dark-energy models 
has two immediate consequences. First, the growth factor is 
changed, which determines how structures grow linearly against 
the expanding background. In models with fixed parameters 
£2o and £2q, structures form earlier if w is larger. Second, un- 
der equal circumstances, the cosmic volume shrinks as w in- 
creases. Dark energy models thus "interpolate" between low- 
density, spatially flat models with cosmological constant and 
low-density, open models. 

For our purposes, we can neglect the clustering of the 
quintessence field and assume that the dark-matter power spec- 
trum is given by the common CDM spectrum. We take the shape 
parameter to be given by F = £2o/i and normalise the spectrum 
such that the COBE-DMR measurements of CMB fluctuations 
on large angular scales are reproduced. This implies a third cos- 
mological consequence. As w increases, the gravitational poten- 
tial of matter fluctuations evolves more rapidly along a given line 
of sight. Secondary anisotropics in the CMB caused by the inte- 
grated Sachs-Wolfe effect thus grow in amplitude. Keeping the 
total fluctuation amplitude fixed to the COBE-DMR data thus re- 
quires the amplitude of the primordial fluctuations to decrease. 
Expressing the power-spectrum normalisation by ag, this im- 
plies that 08 must decrease as w increases. The decrease is gen- 
tle for — 1 < w < —0.6 and steepens as w increases further. 

These findings have two counter-acting effects on the evo- 
lution of dark-matter haloes. First, haloes forming earlier are 
more concentrated because their core density reflects the den- 
sity of the background universe at their formation time. Since 
structures form earlier in dark energy models as w is increased, 
haloes are expected to become more concentrated as w grows. 
Second, the decrease of as with increasing w has the opposite 
effect on the halo formation time and indirectly on halo concen- 
tration. However, cosmologically interesting dark-energy equa- 
tions of state must yield cosmic acceleration today and require 
— 1 < w < —0.6. Within that range, the first effect is dominant, 
and the overall behaviour is monotonic. 

Ideally, extensive, high-resolution numerical simulations 
would be necessary for quantifying the net result of these two 



effects. However, simple algorithms for calculating halo con- 
centrations have already been derived from existing numerical 
simulations. We used them for our work, assuming that they are 
also valid with the modification of Friedmann's equation caused 
by the introduction of dark energy instead of a cosmological 
constant. 

We used three different recipes for computing halo concentra- 
tions. Albeit differing in detail, they agree in concept. Haloes 
are assigned a formation epoch, essentially requiring that a cer- 
tain fraction of the final halo mass has already collapsed into 
sufficiently massive progenitors. The characteristic density of 
the haloes is then taken to be proportional to the mean back- 
ground density of the universe at the halo formation epoch. 
We showed that all three recipes lead to the result that haloes 
are expected to be increasingly more concentrated as w grows 
in quintessence models, showing that the effect of their earlier 
growth is stronger than the effect of decreasing o&. This holds 
for w < —0.6 and reverses for larger w because the integrated 
Sachs-Wolfe effect then requires a steep decrease in c^. The 
particular recipe for computing halo concentrations described 
by Bullock et al. implies that haloes should be ~ 50% more con- 
centrated for w = —0.6 than for w = — 1, where the increase is 
roughly linear with w. 

Finally, we described that halo searches using weak-lensing 
techniques are sensitive to halo concentrations. Using the Sheth- 
Tormen modification of the Press-Schechter mass function for 
quantifying the halo population in mass and redshift, and the 
aperture mass technique for quantifying the weak-lensing ef- 
fects of haloes, we showed that the expected number density 
on the sky of haloes causing 5-a weak-lensing detections ap- 
proximately doubles as w increases from —1 to —0.6, where the 
increase is linear with w . Our results indicate that halo concen- 
trations may be a sensitive probe for the dark-energy equation 
of state, and that gravitational lensing may provide the observa- 
tional tools for applying that probe. 

Note, however, that we did not allow variations in some cos- 
mological parameters which may also change the number of 
weak-lensing haloes. In particular, an effect may arise from 
varying the index n of the dark-matter power spectrum because 
it directly affects the determination of <5%. On the other hand, 
our approach here is to characterise the main effects of dark en- 
ergy on halo formation and to propose weak lensing studies as 
a tool for constraining the dark energy itself, assuming the main 
cosmological parameters will be measured by independent ob- 
servations like those of the CMB. 

Thus, weak lensing turns out to be a powerful tool not only for 
mapping the distribution of matter in the Universe, but also for 
probing fundamental dark-energy properties. Currently, weak 
lensing observations do not allow detailed reconstructions of 
halo density profiles (Mellier 2001, Clowe et al. 2000; Mellier 
& Van Waerbeke 2001), mainly because of the resolution limit 
due to the finite number density of background galaxies. 

On the other hand, interesting new perspectives have been 
opened by several recent wide-field cosmic-shear studies (Ba- 
con, Refregier & Ellis 2000; Bacon, Massey, Refregier, El- 
lis 2002; Wittman et al. 2000, Van Waerbeke et al. 2000, 
Kaiser et al. 2000). Future weak lensing surveys will 
cover even larger fields and, thanks to the improved con- 
trol of systematic errors, they will allow tighter constraints 
on the cosmological parameters. Besides the wide field 
telescopes which are currently in their project study phases, 
we specifically menti on the "dark matter telescope" LSST 



field survey with excellent image quality and depth. 
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( http://dmtelescop.org , propos ed to scan a 7 squa re-degree 
sky field), the VISTA survey (bttp://www. vista.ac.uk ), and the 
SNAP satellite ( frttp://snap. lbl.gov ), whose weak lensing survey 
will cover an area of 300 square degrees, resulting in a very wide 



